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ABSTRACT 

Accretion studies have been focused on the flow around bodies with point mass gravitational potentials, but 
few general results are available for non-point mass distributions. Here, we study the accretion flow onto non- 
divergent, core potentials moving through a background medium. We use Plummer and Hernquist potentials as 
examples to study gas accretion onto star clusters, dwarf and large galaxy halos and galaxy clusters in a variety 
of astrophysical environments. The general conditions required for a core potential to collectively accrete 
large quantities of gas from the external medium are derived using both simulations and analytic results. The 
consequences of large mass accumulation in galaxy nuclei, dwarf galaxies and star clusters are twofold. First, if 
the gas cools effectively star formation can be triggered, generating new stellar members in the system. Second, 
if the collective potential of the system is able to alter the ambient gas properties before the gas is accreted onto 
the individual core members, the augmented mass supply rates could significantly alter the state of the various 
accreting stellar populations and result in an enhanced central black hole accretion luminosity. 

Subject headings: accretion, accretion disks, globular clusters:general, hydrodynamics, galaxies: general 



1. INTRODUCTION 

The process by which stars gravitationally capture ambient 
matter is called accretion. Calculations of the accretion flow 
onto stars are, in general, very difficult. First, one must de- 
termine the flow geometry; in general, if the gas possesses 
intrinsic angular momentum, it will be multidimensional de- 
pending upon the flow symmetry. Second, one must enumer- 
ate the dominant heating and cooling mechanisms that char- 
acterize the accreting plasma. If the gas is optically thick to 
the emitted radiation, the net cooling and heating rates will 
depend on the radiation field. Third, the effect of radiation 
pressure in holding back the flow must be properly accounted 
for. Fourth, one must understand the flow boundary condi- 
tions, both at large distances and at the surface of the accretor. 
Fifth, the possible role of magnetic fields in the plasma must 
be assessed. So, it is not surprising that the problem of gas 
accretion onto stars has been solved in only a few idealized 
cases. In spite of the difficulties, significant progress has been 
achieved in understanding accretion flows by means of such 
idealized solutions. 

The problem of accretion of gas by a star in relative motion 
with respect to the gas was first considered by Hoyle and Lyt- 
deton (1939) and later by Bondi and Hoyle (1944) (BHL ac- 
cretion). The spherically symmetric case, when the accreting 
mass is at rest with respect to the ambient gas, was first studied 
by Bondi (1952) and is commonly referred to as Bondi accre- 
tion. This gives a reasonable approximation to the real situa- 
tion of an isolated star accreting from the interstellar medium, 
provided that the angular momentum, magnetic field strength 
and bulk motion of the interstellar gas with respect to the star 
can be neglected. For other types of accretion flows, such as 
those in close binary systems, spherical symmetry is rarely 
a good approximation. Nonetheless, the spherical accretion 
problem is of great relevance for the theory, as it introduces 
some key concepts which have much broader validity. What 
is more, it is possible to give a fairly exact treatment, allowing 
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US to gain insight into more complicated situations. 

Accretion studies have been focused on the flow around 
bodies with point mass gravitational potentials, but few gen- 
eral results are available for non-point mass distributions. 
Mass models with tractable potentials may accurately approx- 
imate the potentials of galaxies, galaxy clusters and dense 
stellar systems. In this paper we consider calculations of the 
gas flow onto these smooth core potentials. We are interested 
in determining the amount of mass accumulation for core po- 
tentials with a range of velocities relative to the ambient gas. 
In §2, we briefly recapitulate the relevant equations which de- 
scribe the gas structure and deduce analytically the conditions 
for generalized core potentials to accumulate gas collectively. 
Because the gas structure is multi-dimensional in nature, we 
adopt a numerical approach to verify our analytic solutions. 
We describe the numerical scheme and the range of model 
parameters in §3, and the results of these calculations in §4. 
These results are applied, in §5, to study the accumulation of 
gas into dense stellar systems, galaxies and galaxy clusters 
moving in a variety of environments. In §6, we summarize 
our results and discuss their implications. 



2. MASS ACCUMULATION IN A CLUSTER CORE POTENTIAL 
2. L Generalized Solutions 

Much of the effort herein will be dedicated to understand- 
ing the conditions by which smooth core potentials are able 
to accrete surrounding gas with highly enhanced rates. In 
particular, we would like to understand under which condi- 
tions a stable density enhancement will persist within a given 
potential well and derive a relation between the gas relative 
velocity, the local sound speed and the the central velocity 
dispersion of the cluster potential, which can be carried over 
quite generally to more complicated accretion flows. 

To treat t he pro blem mathematically we proceed as in |Lin| 
& Murray ( 2007j i. We take spherical polar coordinates with 
origin at the center of the potential. In spherical symmetry, the 
fluid variables are independent of angle and the gas velocity u 
has only a radial component. For a steady flow, the continuity 
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equation reduces to 



r pu= constant. 



(1) 



This constant is commonly denoted to as M, and refers to a 
mass inflow (outflow) rate. In the Euler equation the only 
contribution to the external force comes from gravity and this 
has only a radial component 



du 1 dP dd> 
dr p dr dr 



(2) 



Hence, the term {l/p)i^) is (cl/p){dp/dr), where c, is the 
sound speed. We replace the energy equation by the poly- 
ti'opic relation 



P = Kp^, constant. 



(3) 



In this case: cl = dP/dp = K'^p''~^ . This allows us to approx- 
imately treat both adiabatic (7 = 5/3) and isothermal (7 = 1) 
gas flows. After the solution has been found, the adiabatic or 
isothermal assumption should be justified by consideration of 
the particular radiative cooling and heating of the gas. 

The gravitational potential of dense stellar systems, galax- 
ies and galaxy clusters can be accurately described by smooth 
core potentials with 



4>(r) = - 



90 



(4) 



where (Sc is the radius of the core. The potential's core scale 
is given by 

C(r)=(r/' + a3f)" (5) 

Cores potentials described by Q are non-divergent, so that 
(p{r > 0) = constant. This property, as we will show, allows 
core potentials to sustain subsonic inflow and collect ambient 
gas into a quasistatic envelope. 

An estimate of the radial density distribution of collected 
ambient gas, p(r), follows from the application of the conser- 
vation laws for energy and momentum. From equation Q we 
have 



dp u cl 
p 27-1 



+ ffl = const. 



(6) 



Equation (|6]l makes use of the polytropic equation of state, 
where dp = K^p^'^dp and c\ = ^p j p = K^p^'"^ . Here, the 
limits (f){r i-> 00) = and u{r ^ 00) = help determine the 
constant 

+ 0=^. (7) 
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Dividing equation (j7| by cl and making use of c^/c; 



equation in Bondi ( 1952 1 



/p2o', one recovers a modified form of the Bernoulli 
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(8) 



The isothermal version of the Bernoulli equation can be de- 
rived from equation (|6| using dp = Kdp and c^ = c^, 



1m" , , 

-^ + logp+ — =logPo 



(9) 



Together with knowledge of the radial velocity profile, 
equation ([8]l allows one to determine p(r). For subsonic in- 
flow u^(y)/cl ^ 1, the ambient gas collects into a a quasistatic 
envelope with a density distribution given by 



p{r): 



1 + 



(7-1)00 



(10) 



where poo is the density of the ambient medium. For an 
isothermal equation of state, this result is simply 



p{r) = Poo exp 



(11) 



In 5 2. 1 . 1 we derive under what conditions we can safely 
ignore the uirf' /cl term in equation ( 10 1. 



2.1.1. Subsonic Flows in Core Potentials 

We now consider the validity of the subsonic inflow as- 
sumption \u(r)\ < Cs used to derive equation (fTo| assuming 
subsonic flow at infinity, |moo| < c^o- In spherically symme- 
try, the radial structure of the flow's velocity can be calculated 
by integrating (|2]i with the help of ([3]) and ([T]i, so that 

At large distances from the potential the term (2cl-r^) on 



the right hand side of equation (Il2b must be positive since cl 
approaches and cfxx 1/r , while r increases without limit. 
Given the that the inflow is initially subsonic, it will remain 
subsonic provided 

2 r d<j) 



ct > 



2 dr 



(13) 



for all r in the flow such that there is no transonic point. 
Because Cs is a monotonically decreasing function of r, 



equation (13 1 is satisfied provided 



9 r ua> 

c > 

°" 2dr' 



(14) 



Consequently, condition (14 1 guarantees \u(r)\ < c^ through- 
out the gas flow. 



Describing 0(r) as in equation (|4|i, condition ( 14 1 can be 
rewritten as 

1 /3ar^0o 



cL > 



The right hand side of equation ( [T5| has a maximum at 



(15) 
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SO that condition ( 14 1 is satisfied for core potentials if 

00 /3 



(— 

\a+l 



(16) 



(17) 



Notice that as U5c 0, the ambient gas is unable to collect 
into a quasistatic envelope and, as a result, the density struc- 
ture of the flow is no longer accurately described by equation 



10 1. In this limit, the potential essentially reduces to that of a 



pomt mass and the inflow rate can be derived following Bondi 
( )1952) . 
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2.1.2. Mass Accumulation in Non-Stationary Core Potentials 

The density structure of ambient gas collecting into a core 
potential in the case that the gravitational structure is col- 
lecting gas isotropically is described by equation ( [TO) l. If 
the potential is moving with some velocity with respect to 
the ambient medium, one expects the build up of a stable 
density enhancement to be inhibited. A characterization of 
this impedance on the growth of the central density enhance- 
ment is derived here by comparing the pressure support due 
to the collected ambient gas to that provided by the ram pres- 
sure. The radius at which these two pressures are equal de- 
fines a characteristic radius, which we denote here as r^. 
A schematic diagram illustrating the geometry of the flow 
around non-stationary core potentials, together with a sketch 
of the expected sound speed and density profiles of the col- 
lected ambient material, is given in Figure T] 

The characteristic radius, r^, is calculated here under the 
assumption that a stable density enhancement has been es- 
tablished, whose radial profile is well described by equa- 
The potential, together with its collected ambi- 



lOi 



tion 

ent gas, is then assumed to move with a constant velocity, 
"const ~ Moo = Moo Coo- This providcs a reasonable approxima- 
tion provided that the flow variables across the shock region 
are accurately described in terms of their limiting values (i.e. 
the incoming flow is not significantly altered by the gravita- 
tional potential and therefore the interior flow preserves its 
static solution). When this high velocity gas hits the station- 
ary ambient medium, a bow shock forms. In the frame moving 
with the shock, the density jump condition can be written as 



(7+1)/^; 



(7+l) + (7-l)(AiL-l) 



(18) 



where = p{r), with p{r) given by equation JTOb. Equation 
( [T8] l assumes the adiabatic index 7 is the sarneror the post 
and pre-shock gas. The jump conditions across the shocked 
region yield 



(/io(7-l) 



(7+l¥5o 
2 + (7-lWo 



7-1 



-1 



(19) 



where we have used Q{rs) = (rf + 03^)". These standoff dis- 
tance define the characteristic radius of the bow shock. For 
00 A# we have a- 13= 1, so that 



C(r, = OJc) = 2 



^91/" 



(20) 



Equation (19 1 can be rewritten for an isothermal equation 
of state as 



00 1 



for 7^1. 



(21) 



In general, core potentials moving at high velocities have 
density enhancements and characteristic radii that decrease 
with increasing velocity. At sufficiently high Mach numbers, 
one expects the incoming flow's ram pressure to overcome 
the central pressure due to the collected ambient gas and to 
penetrate into the potential's core. The transition between a 
stable density enhancement and an unsupported, supersonic 
inflow occurs roughly when w CSc. When < XS^ , the abil- 
ity of the core to sustain a stable density enhancement is lost. 
For this reason, we expect rapidly moving potentials to have 
smaller density enhancements than those predicted by equa- 
tion ( [T0| . In these cases our formahsm inherently overesti- 



mates the internal pressure support against the inflowing am- 
bient material, which in turn overestimates the characteristic 
radius. The detailed behavior of the characteristic radius with 
adiabatic index, Mach number and gravitational potential type 
will be explored in fusing numerical simulations. In the 
following numerical models, we will approximate isothermal 
equations of state by setting 7 = 1.01. For consistent compar- 
ison, we will use equations ( 10 1 and ( 20 1 with 7 = 1.01 in all 



plots unless otherwise noted. 

2.2. Specific Core Potentials 

In this paper we focus our attention on the ability of stellar 
clusters, galaxies, and galaxy clusters to collect ambient mate- 
rial. Stella r clusters can be ac c urately described by Plummer 
poten tials ( Briins et al. 2009 [ [Pflamm-Altenburg & Kroupa 
|2009l l 



0(r): 



GMp 

(r2 + r2)i/2' 



(22) 



where Mp is the mass of the potential and Tc is the commonly 
referred to as core radius. In this case, 0o = GMp, XS^ = r^, 
(3 = 2, and a = 1/2, so that equation ( 10 1 can be rewritten as 



P(r) = Po 



1 + 



GMp(7- 1) 



(r2 + r2) 



1/2 



(23) 



Dark matter profiles of galaxies and gal actic clusters are, 
on the other hand, commonly modeled by Hernquist ( 1990| l 
potentials 

0(.)=^, (24) 



r+a 



where Mj, is the mass of the potential and a is the scale length. 
In these cases, 00 = GMh, XSc = a, f3 = 1, and a = 1, so that 
equation ( [T0| becomes 



P(r) = Po 



1 + 



GMh(7-l) 



(25) 



In what follows, these two types of potentials are numerically 
modeled moving through ambient gas at rest to and compare 
to our analytic treatment. 

3. SIMULATIONS OF MASS ACCUMULATION IN CORE 
POTENTIALS 

3.1. Numerical Method and Initial Model 

To examine the ability of core potentials to collect ambient 
gas into a quasistatic envelope, we simulate core potentials 
moving through ambient gas with FLASH, a parallel, adaptive 
mesh refinement hydrodyna mics code. This sche me, and tests 
of the code are described in Fryxell et al. ( 2000| l. 

Common to all calculations is the placement of inflow 
boundaries to one side of our computational domain in or- 
der to simulate the motion of the core potential through the 
external medium. As a result, ambient gas moves to the right 
of our grid at a supersonic velocity. When this high velocity 
gas hits the stationary potential, a bow shock forms. All grav- 
itational potentials are modeled here using a Plummer ( [22| or 
a Hernquist ( [24| functional form. 

Simulations start with a uniform background density and 
run until a steady density enhancement is established inside 
the core's pote ntial. This usually t akes about 10-100 sound 
crossing times ( Naiman et al.| l2009). Several models were 
run longer to test convergence and density enhancements were 
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found to change only slightly with longer run times. We fur- 
ther tested convergence of our models for several resolutions 
and domain sizes. All tests produced similar density enhance- 
ments to those shown here. After hundreds of sound cross- 
ing times, the flow is relatively stable, and does not exhibit 
the flip-flop instability seen in two dimensional simulations 
(jBlondin & Pope 2009). 

The effects of self gravity of the gas are ignored. This is 
adequate for most of our models, for which the accumulated 
mass is significantly less than the mass responsible for the 
potential. To improve the controlled nature of the models, and 
to facilitate comparison with our analytic estimates, we do not 
explicitly include radiative heating or cooling. As with our 
analytic estimates, we assume a polytropic equation of state 
P (X p'' . The effects of radiative equilibrium are approximated 
by having the gas evolve with an adiabatic constant 7 = 1.01 
(giving nearly isothermal behavior), while inefficient cooling 
is model here by having the gas evolve with 7 = 5/3. 

3.2. General Properties of the Gas Flow 

The classical Bondi treatment for non-divergent potentials, 
representative of star cluster and galactic gravitational poten- 
tials, is only a fair approximation far from the core when 
GMjm^ > +M^. When GMjUi^ ^ the col- 

lective potential alters the local gas properties, allowing the 
core to sustain subsonic mass inflow and, as a result, collect 
ambient gas into a quasistatic envelope. 

Figure |2] shows the resulting density structures of the col- 
lected ambient gas for a point mass and two different core po- 
tentials. To facilitate comparison, the Plummer and Hernquist 
potentials depicted in Figure [2 have been constructed to yield 
the same circular velocity pe^ radius: r^. This requirement 
can be rewritten as 

a = Vlr^ , (26) 

which, by demanding the mass within this characteristic radial 
scale to be the same, gives 



Hence, equations ( 23 1 and ( 25 1 give 



Ph(r = 0) ( 8 



Pv(r = 0) V3%/3 



(27) 



(28) 



where p^^^ = p{r = 0) > Poo- 

The flow pattern around a core potential is multi- 
dimensional and complex. In the frame of the potential, the 
gas streamlines are bent towards the cluster center Some shall 
intersect the center, while others converge along a line behind 
it. The convergence speed of the gas determines the reduc- 
tion in its velocity relative to the potential due to shocks, and 
therefore whether or not the gas is accumulated in the core. As 
seen in Figure|2] the flow structure within the cluster markedly 
differ from the classical Bondi flow. 

The properties of the core potential have a decisive effect 
on the radial profiles of the collected ambient gas. For r ^ CJc 
and Pmax Poo. equation (flOjl gives 



pjr) 

Pmax 



\-a{rl^,y 



(29) 



From equations (29i and (26 1, it follows that p/pmax oc [1- 

(r/\/2r(,)]~ and p/ Pmio. [l-2(r/rc)^]^ for Hernquist and 
Plummer potentials, respectively. Figure [2] shows that the the 
steeper Hernquist interior potential produces gas accumula- 
tion which is more centrally confined within its core than that 
in the Plummer potential. 

At r ^ n3c, the divergence of the gravitational potential nat- 
urally leads to a loss of pressure support against the incom- 
ing material. As a result, the characteristic radius, rj, will 
be smaller for less extended (or less massive) core potentials. 
Inserting equation (127]) into equation ( 19 1 and making use of 
equation ( 20 1, one finds 



C(''s = roc)h 8 



CC'-s = roc)p 3V3 



(30) 



This shows that gas collected in potentials with steep core 
profiles have, as expected, sharper central density profiles. 



for Hernquist and Plummer potentials yielding the same cir- 
cular velocity peak radius, 7, /ioo and Cqo- Under these condi- 
tions, the ambient gas collected into a Hernquist potential at 
r':^Uic will be comparatively denser and thus lead to an in- 
crease in pressure support against the incoming flow. This 
is mainly due to the comparatively larger mass (see equa- 
tions ( [27| ) and ( [T9| ) together with the less divergent poten- 
tial at large radii. This enhanced support forces the flow to 
move around the central density enhancement, thus leading to 
a widening of the standing bow shock. This is evident in Fig- 
ure |3] which plo ts the bow shock patterns produced around a 
point mass (Ruf fert|1994 1, Plummer and Hernquist potential 
set in motion through an initially uniform medium. 

3.3. Mach Number Dependence 

In line with the conventional treatment, potential cores 
moving with respect to the ambient medium at increasing su- 
personic velocities will have density enhancements that are 
progressively smaller and significantly more offset from the 
cluster's center. Figure |4] shows the resulting density and ve- 
locity flow structures generated around a point mass, Plum- 
mer and Hernquist potential moving through an initially uni- 
form medium with /ioo > 1 ■ The differences in the flow struc- 
ture are most notable when comparing the flow velocity struc- 
tures depicted in Figure]?] which show that the central density 
enhancements in core potentials are quasi-static, while those 
for the point mass are not. While there are regions of in- 
flowing gas around core potentials, primarily behind the bow 
shock and downstream from the core's center, there are large 
central regions that experience no net inward gas flow. 

A stable density enhancement is observed to persist as long 
as /ioo ^ 3 (Figure]?]). For larger relative velocities, the central 
gas pressure support is overwhelmed by the ram pressure of 
the incoming gas and, as a result, the bow shock begins to 
detach from the potential's core. The transition from a bow 
to a tail shock takes place at around /ioo ~ 3. Note that in 
the case of a point mass potential with an arbitrarily small 
sink radius, the gravitational forces will always overwhelm 
the ram pressure at a fini te radius. In th is case, there is no 
transition to a tail shock ( Wolfson|1977| l. However, for core 
potentials (and point mass potentials with a finite sink size), 
the gravitational force approaches a constant value near the 
potential's center, thus allowing the formation of a tail shock. 

As shown in Figure ]5] the radial density profiles of the ac- 
cumulated ambient gas, within a core potential differ from the 
classical Bondi profiles, and, for low Mach numbers, are ac- 
curately described by the quasistatic solutions derived in ^ 
These solutions fail to provide a good description of the den- 
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sity of the inflowing gas for ^ 1- This is because these 
solutions were derived assuming a static ambient gas. The de- 
crease of the characteristic radius and of the collected ambient 
gas density with increasing /ioo are clearly visible in the sim- 
ulations. This decrease is connected with an increase of the 
ram pressure of the incoming material which, in turn, helps 
overcome the central pressure of the collected ambient gas. 

3.4. Dependence on the Polytropic Index 7 

The determination of the internal structure of gas flowing 
around gravitational potentials depends upon detailed knowl- 
edge of the equation of state. In ^ we have derived analytical 
solutions for the internal pressure and density structure of the 
flow in the cases in which the pressure depends solely upon 
the density through an equation of state of the form P = Kp'^ . 

The dependence of the overall gas density structure on the 
polytropic index 7 is illustrated in Figures |6 and|7] By com- 
paring the flow structures in Figures |6] and 7] one sees that 
ambient gas is in fact collected much more gradually when 
its pressure support increases adiabatically (7 = 5/3). From 

I I — r 

equation ( 10 1, it follows that P cx ^ . Thus, gas accumulated 
in potentials with higher adiabatic indexes will be less cen- 
trally confined and, as a result, will be less effective at protect- 
ing the collected gas from being displaced by the inflowing 
material. By contrast, gas collected isothermally (7^1) into 
a core potential will be comparatively denser and thus lead to 
a significant increase in pressure support at small radii. This 
is consistent with the behavior observed in Figure l6j where 
decreasing 7 moves the shock closer to the potential's core. 
Despite the initially less extended bow shock, the enhanced 
support provided by more compressible flows at small radii 
delays the transition to a tail shock as the relative velocity in- 
creases and, as a result, the bow shock radius decreases. 

3.5. General Conditions Required for Large Density 
Enhancements 

A core potential can sustain a large central density enhance- 
ment provided that rs > tDc. This condition can be rewritten as 
GM/n3c > (c^ + M^) = y^, or alternatively, as Cv V where 
(Tv is the velocity dispersion of the system. 

In Figure |8] we examine the ability of a core potential to 
collect ambient gas as a function of the characteristic radius 
Tj. A density enhancement is observed to persist as long as 
Ts > 03c. Once < tDc, the collective potential is unable to 
alter the properties of the flow within r < CJc . 

In addition, the equation of state alters the compressional 
properties of the incoming flow. Isothermal flows, corre- 
sponding to 7 « 1, allow for the effective build up of mass. 
Adiabatic gas, 7 = 5/3, will instead rapidly heat the gas as 
it accumulates in the core, and halt the gas inflow at much 
lower central densities than for the 7 w 1 case. Figurel9]com- 
pares the simulated density enhancements results with those 
predicted by equation ([TO]l. 

3.6. Bow Shock Instabilities 

Most shock instabilities identified and studied in astro- 
physics are related to the acceleration or deceleration of the 
shock itself (e.g. Bernstein & Book 1978, Vishniac 1983). 
By definition, stationary shocks are stable with respect to 
these mechanisms, but can still be unstable. Nakayama (1993) 
pointed out the radial instability of the shock if the flow is im- 
mediately accelerated after the shock transition, in isothermal 



flows. In this respect the shocked Bondi flow should be un- 
stable. The validity of the postshock acceleration criterion for 
adiabatic flows is still uncertain, even for radial perturbations 
(Nakayama 1994). Independently of the sign of post-shock 
acceleration, Foglizzo & Tagger (2000) and Foglizzo (2001) 
revealed another generic instability mechanism based on the 
cycle of acoustic and vortical perturbations in the shock re- 
gion. 

The vortical-acoustic instability is fundamentally non- 
radial. In this case, both vorticity and entropy perturbations 
are advected from the shock to the accretor, and both are cou- 
pled to the acoustic perturbations. The most unstable cycle 
involves high frequency acoustic waves, those able to explore 
the hottest parts of the flow but still be refracted out, with a 
wavelength slightly larger than the smallest size of the sonic 
surface. In the isothermal limit (7 1), where entropy pertur- 
bations are absent from the problem, vorticity is more appro- 
priate than entropy to describe the advective-acoustic cycle. 

The entropic-acoustic cycle is expected to be an efficient 
instability mechanism for 7 « 5/3, provided /loo ?J 2 and that 
the distance to the sonic surface (with depends on r^) is small 
enough (Foglizzo et al. 2005). This coincides with our numer- 
ical simulations of core accretion which show strong instabil- 
ity when ss Ui^ (Figure [8]l. The lack of strong instability 
of adiabatic flows around core potentials with rs ^ CSc could 
be due to the lack of effective acoustic feedback from the ad- 
vected vorticity perturbation. Our numerical experiments also 
show the presence of instabilities within fast moving, isother- 
mal flows. In particular isothermal flows appears much more 
unstable when r^ « XS^ (Figure l6]l. In this conditions, a small 
pressure perturbation of the shock is able to generate vorticity 
perturbations very efficiently for 7=1 strong shocks. The dif- 
ferences seen in the shock geometry between these two lim- 
iting cases precludes an accurate comparison of the stability 
of the flow, especially since the instability threshold is very 
sensitive to the size and shape of the sonic surface (Foglizzo 
et al. 2005). 

4. ASTROPHYSICAL RELEVANCE 

The ability of core potentials to collect ambient mass de- 
pends sensitively on their velocity dispersion, relative speed 
with respect to the ambient gas and gas cooling properties. 
Estimates of velocity dispersion and environmental gas prop- 
erties are known for a wide range of systems, including star 
clusters, dwarf galaxies, large galactic halos, and clusters of 
galaxies. In what follows, we refine our calculations to in- 
clude a wide variety of embedding media and discuss how 
environmental properties conspire to enable or disable the ef- 
fective collection of surrounding mass in these systems. We 
assume the individual stellar members are accreting at low 
enough rates such that mass removal (and feedback) can be 
neglected. These results are summarized in Table [T] 

4. 1 . Star Clusters 

Star clusters are gene rally well de scribed by Plummer po- 
tentials (Plummer \9\\ \ Briins et al. 2009; Pflamm-Altenburg] 
'& Kroupa 2009). In this section, we will discuss two lim- 
iting cases spanning the range of environments star clusters 
might inhabit. First, we discuss a compact star cluster mov- 
ing slowly with respect to cold gas, representative of young 
clusters in a galaxy merger environment. This situation is sim- 
ulated in Figure 10 (model A). Second, we model a less com- 
pact cluster moving quickly with respect to the background 
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gas. These less favorable conditions for gas accumulation are 
simulated in Figure 10 (model B). 



4.1.1. Young Star Clusters 

Young massive star clusters are frequently found in merging 
systems lik e the Antennae and Ca rtwheel galaxies (Whitmore] 
S al. 1999; |Crivellari et al.| 2009|l. In the Antennae complex, 
for example, young (compact) star clusters are observed to 
move slowly w ith respect to the cold ga s: « 10km s"' 
and /j,oo ^ 2 (!i |Gilbert & Graham 2007 1. Observations in 
the infr ared suggest large amounts of dust present in the gas 
(Brandl et al. 2005 ), which justifies an assumption of isother- 
mal gas. A typical cluster, Mp = 3.5 x 10^ M©, OJc = /"c = Ipc, 
in this cold, efficiently cooling gas (7 = 1.01), moving at 
/ioo = 2.0 is shown in Figure [10| (model A). In these merg- 
ing environments, conditions are ideal for accumulating large 
amounts of ambient gas in these systems as ^ r^.. Tak- 
ing Poo ss 10~^^g cm"', similar to that inferred by Zhu et al. 



4.2. Galactic Systems 

Hernquist potentials (| Hemquist|1990 |l are used to describe 
a wide variety of galactic systems - from the dark matter ha- 
los of dwarf galaxies an d Milky Way sized systems, to the 
halos of galaxy clust ers ( [Aden et al.||2009[ [Scott et al.||2009 



pii 



[Hwang & Lee|[2008] l. In what follows we assume, for sim 



pTicity, that these gravitational potentials are not significantly 
modified from their dark matter profiles by the presence of 
gas or stars. 

4.2. 1 . Dwarf Galaxy Halos 

Dwarf galaxies (DGs) are seen eithe r moving through the 
IGM or embedded within larger halos ( Diemand et al.|2008 



Mateo|1998 ). For those moving under the influence of a larger 



( |2003[ l, we find the collected mass to be a sizable fraction of 
the total cluster mass: M^^^/M^ k. 0.05. 

For these young clusters (< Gyr), one might worry about 
the applicability of our results, as feedback due to massive 
stellar winds is neglected in our treatment. We can estimate 
the magnitude of this effect by comparing the ram pressure 
from stellar winds of the cluster members to the pressure sup- 
port generated by the collected ambient gas. The ram pres- 
sure due to expanding stellar winds is approximately given by 
P* = /9* = NM^, \\ I (Anr^) for massive stars. Each star here 
is assumed to have a mass loss rate and wind velocity and 
V*, respectively, expanding into a spherical volume of radius 
r. We can then calculate the size of the wind bubble, r^, by 

equating P*(rw) = Pp(r„), where Pp = c^[p(r)]'^/7 R:!c^(p), ^ in • r- 

where (p) is the average density in the cluster's core. Solv- lor model D m Figure 



ing for r„ gives « 0.02A^2 ^ 



,1/2 



-1/2 



■18 where 



gravitational potential we expect that, on average, DGs to or- 
bit around the ambient gas with a velocity comparable to the 
circular velocity of their host: > GMh/a. However, many 
DGs are found to reside in nearly radial orbits, and thus will 
spend only a small fraction of their orbital time moving at 
lower velocities (Die mand et al . 2008 2007). 

For example, a DG on a highly elliptical orbit, e = 0.9, will 
spend approximately 1 % of its orbital period moving at less 
than 20% of the host's circular velocity. During this period a 
large central density enhancement is expected to be effectively 
sustained within the DG' core ifGMh/fl>c^+M^. That is, if 
the DG is surrounded by cold gas with c^^ < GMh/a. Such an 
ideal situation for mass accumulation is illustrated by model 
C in Figure pT| In these conditions, if the density of cold gas 
is Poo = 10 gem"-"*, the mass collected by the DG could be as 
large as^Macc/Mh ~ 0.14. In contrast, if the surrounding gas 
is virianzed, Cqo ~ cTv.hoiit ^ GMh/a, the DG halo will spend 
most of its orbit with mass inflows similar to those depicted 
In this case, for poo = 10"^^ gem"-', 
the accumulated mass is only M^^^/Mh « 5 x 10"^. 



= 102a^2, M* = 1O"^M*,-7M0 yr"i, = 10^v*,3 kms"', 
(p) = 10"'^ (p)-i8gcm"-' and Cqo = lOcoo.i km s"'. For large 
density enhancements, we have ^ /"c « 1 pc and, therefore 
the stellar winds would be unable to significantly impede the 
inflowing of gas. 

4.1.2. Globular Clusters 

Globular cluster' s orbits are primaril y dictated by the host 
galactic potential ( |Harris| [19951 |1991| l. In the Milky Way, 
globular clusters (GCs) are obser vea to move at highly su- 
personic velocities poo ~ 5 - 10 (jDinescu et al. 1997| l. In 
addition, the gas within the galactic halo is expected to 
be virialized so that Coo ~ Uoo, implying ^ GMp/CSc. 
As a result, for a typical GC with Mp = 3.5 x lO^M© and 
tDc = '"c = 2pc moving at poo = 4.0, one expects negligible 
mass accumulation as < r^. This is true even if the clus- 
ter is surrounded by cold, efficiently cooling ambient gas 

(Co 



4.2.2. Halos in Merging Systems 

A massive galactic halo is likely to experience a relatively 
equal mass merger during its lifetime. Galaxies of M w 
10" M0 will have undergone, on average, about one merger 



with a galaxy of M > 5 x 10 M0 since z = 1.2 (Robaina et 
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(model 

B). In this case, for a background density similar to that of 



lOkms and 7 = 1.01) as shown in Figure 



ar]|2010 ). Because the smaller component is moving under 
me gravitational influence of a halo of similar mass, we ex- 
pect Moo ~ o'v Moreover, if gas is effectively virialized within 
more massive halo, then Coo ~ <Jy and, as a result, a large 
density build up is expected to be effectively retained in the 
core of the smaller halo. This situation is depicted in Fig- 
ure 12 (model E), in which a a Milky Way-like sized halo 
is simulated moving through hot gas with Cqo = lOOkm/s 
at Poo = 2.0. Comparing the bremsstrahlung cooling and 
the core sound crossing time, we find fcooi ^ fdyn- Since 
cooling is expected to be inefficient, we assume 7 = 5/3. 
While r, > (53^ = a, ineffective cooling results in an accumu- 
lated mass of only M-^^c/Mh = 10"^ for ambient densities of 
Poo = 10""' 



'gem 



flie ISM, 



10 ^"^gcm the accumulated mass is only 



Macc/Mp « 2 X 10"^ 

A GC might collect gas more effectively if its orbit hap- 
pened to lie in the disc plane of its host galaxy. If the GC mov- 
ing with Poo < 2 through the galactic disc, we might expect 
the mass accumulation to more closely resemble that seen for 
mode l A in Figu re [TQ| However, these conditions are rare for 
GCs ( [Dinescuet ai:.T997) . 



4.2.3. Galaxy Clusters 



Similar to the case discussed above, a large central gas den- 
sity enhancement is expected in equal mass merging galaxy 

^ We note here that model C has been run until a quasi-static density profile 
has been established, which takes about a few core sound crossing times. If 
the DG is at apogee for a shorter time, the accumulated mass is expected to 
be lower. 
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clusters, such as the Bullet cluster complex. The Bullet cluster 
is a gala xy cluster merging with another slightly more massive 
cluster ( Markevitch et al.|2002 l. In this case, Mh = lO'^M©, 
!!!(; = fl = 1 Mpc, fioo = 2.5, and Coo = 1580kms"' . For an inef- 
ficiently cooled gas (7 = 5/3), this results in « 03c = Such 
situation is depicted in Figure [13] (model F). Assuming IGM 
densities of p^o = 10~^^ gcm"^ results in an accumulated mass 
of M.^cjMii = 3.5 X 10"'*. This simple model of the Bullet 
cluster gives similar flows structures to those found in more 
sophisticated calculations which include the stellar and gas 
mass distributions. 

5. DISCUSSION 

In ^we have shown the conditions required to collectively 
accumulate large quantities of gas from the ambient medium 
in the centers of a variety of astrophysical systems thought to 
be well described by a core gravitational potential. The conse- 
quences of large mass accumulation in such systems fall into 
two basic categories. First, if the gas cools effectively star 
formation can be triggered, generating new stellar members 
in the system. Second, if the gas accretes efficiently onto the 
core's current stellar population, the state of these stellar sys- 
tems (and possibly the central black hole) can be significantly 
perturbed. In what follows we discuss the implications of 
gas accumulation and their accompanying observational sig- 
natures. 

5.1. The State of the Accumulated Gas 
5.1.1. Emission Properties 

The emission properties of the accumulated gas in clusters 
of galaxies and young stellar clusters are briefly reviewed in 
this section. At temperatures > 3 x 10^ K, the primary cool- 
ing process for intracluster gas is thermal bremsstrahlung, 
while for temperatures < 3 x 10^ K, line cooling becomes 
very important. For the the intracluster gas in most star clus- 
ters, the cooling time is longer than the Hubble time (age 
of the universe). Thus cooling is not dynamically very im- 
portant in these cases. Figure [14] (bottom panel) shows the 
surface brightness profile, calculated assuming optically thin 
cooling resulting from bremsstrahlung and collisionally ex- 
cited metal transitions, for a cluster merging environment 
with parameters selected to mimic the Bullet Cluster com- 
plex (model F). We find an integrated bolometric luminos- 
ity of Lboi ~ 10'*'* ergs"' and a bowshock structure that are 
comparable to those observed by Chandra (? Markevitch|2006 



[Markevitch et al. 2002). In the center of the cluster, where 
cooling is important, our calculated emission is considerably 
harder than found in s imulations including the c ooling term 
in the energy equation ( jMilosavljevic et al.|2007| l. 

Also shown in Figure |14| is the surface brightness profile 
of the accumulated gas m a young star cluster with param- 
eters similar to those seen for the stellar cluster members in 
the Antennae galaxy merging system (model A). To calcu- 
late the emission we assume optically thin cooling, where 
we employ an effective cooling curve. In the cluster's core, 
the near infrared (NIR) surface brightness is approximately 
« lO'*' ergs"' pc"^, which is probably optimistic given that 
we do not explicitly include the cooling term in the energy 
equation but instead use an adiabatic simulation. This simple 
model gives, nonetheless, results that are consistent with the 
observed range of lumin osities in the Anten nae star clusters: 



Lnir = lO^"*- 10^*2 ergs"' (|Brandl et al. 2005|). For clusters em- 



bedded in these types of mergmg environments, the emission ^^^^^ orders of magnitude. 



from the bowshock and cluster core should distinguishable^ 
given enough spatial resolution and provided the reprocessed 
stellar radiation is not the dominant contribution to the NIR. 

5.1.2. Enhanced Central Densities Leading to Star Formation 

If the gas collected by these potentials cools efficiently, 
we would expect it to fragment and form stars if the den- 
sity is high enough. In general, one can relate the gas den- 
sity to the star formation rate (Kennicutt 1998): Esfr ~ lO"'* 
(Sgas/lO'*M0pc"^)' '* Moyr"'pc"^, where Esfr and Sgas are 
the star formation and gas column densities, respectively. A 
dwarf galaxy, for example, embedded within a larger halo 
with properties similar to those assumed for model C (i.e.. 
Poo = 10"^'*gcm"-' and a 200pc star forming region) is ex- 
pected to form stars at a rate of Estv = 10"' Moyr"', similar 
to tha t predicted by star formation simulations in these sys- 
tems ( |Sawala et al.|2010| ). 

For the star clusters embedded in the Antennae galaxy 
merging system, whose conditions are relatively well de- 
scribed by our model A, we find Ssfr ~ 4 x 10"^ yr"' per 
system, or a total rate of S^fr « 4(A/c/lO^)M0yr"' for all A^c 
clusters. Such star formation rates were calculated assuming a 
star forming region of 5 pc and a background ambient density 
of Poo = 10"^' gem""*, which results in an average increase in 
cluster gas core density of (p) = 10"'^gcm"^ (model A). The 
resulting rates are similar to those observed in the Antennae 
complex ( [Gilbert & Graham||2007| l. If we assume that such 
rates are effectively sustained over the calcula ted orbital pas- 
sage time scale of about 20 Myrs ( Karl et al.^20 10), one ex- 
pects at least 10^ of new stars to be created per system. 
More detailed simulations that take into account self-gravity, 
realistic cooling and feedback effects are needed before de- 
tailed comparisons can be made with the inferred conditions 
required to explain subsequent episodes of star formati on in 
stellar c lusters (Conroy & Spergel 2010, Gratton et al. j2004[ 
[Carretta et al.,2010,,Piotto 2009;Mackey et al..2008) . 

5.2. Accretion Onto Individual Stellar Members and the 
Central Massive Blackhole 

A compact star (or a black hole) of mass M*, moving 
with relative velocity u through a gas of ambient density p^o 
and sound speed Cqo, nominally accretes at the Bondi-Hoyle- 

Lyttleton rate: M 47r(GM*)^poo(M^ + c^)"''^^- If, however, 
the collective potential of the system (e.g. star cluster or 
galaxy) was able to significantly increase the surrounding gas 
density prior to being accreted onto the individual neutron star 
members, then the accretion rates of the individual core stel- 
lar members (or the central black hole) would be enhanced by 
^ (p) I Poo, where (p) is the average increase in gas density in 
the core. It is to this problem that we now turn our attention. 
In particular, we study the effects of enhanced accretion rates 
onto the white dwarf core members and the central massive 
black hole. 

5.2.1. White Dwarfs in Stars Clusters and Galaxies 

In mature stellar systems, we expect a significant fraction 
of the stellar mass to be in white dwarfs. If the collective 
potential of the system (which contains the white dwarfs) is 
able to alter the ambient gas properties before the gas is ac- 
creted onto the individual stars, then the resulting enhanced 

^ This is because the change in brightness across these two components is 
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accretion rates could significantly alter the state of the accret- 
ing stellar population. If, however, the density is not signif- 
icantly increased, then the stars accrete gas as though they 
move through the external medium independently. 

The resulting mass accretion rates for the white dwarf mem- 
bers depends on the radial distribution of both stars and gas. 
Here we calculate the expected mass accretion rates in the 
core of slowly moving star cluster (model A) using two ex- 
treme examples for the radial distribution of white dwarfs 
remnants. The first one is based on Fokker-Planck mod els of 
a core collapsed (centrally condensed) globular cluster ( |Dull| 

tal. 1997), and the second simply assumes that the white 
varfs, containing 1 % of the total mass, follow the radial stel- 
lar mass distribution. FigurefTSlshows the mass accretion rate 
distribution of the white dwarfcluster members for both radial 
distributions (solid and dashed lines respectively). 

This range of mass accreti on rates are expected to re- 
sult in enhanced novae rates ( |Nomoto et al.|[2007| l. With 
a nova outburst rate of '^ra ~ lO"'* yr"\ we expect about 
« O.U^'nr/lO-'^ yr"') (A^„d/10'*)(A^c/102) outbursts per year. 
This is probably a conservative estimate as the individual no- 
vae rates are expected to increase to ^ < 10"' yr~' for th e 
larger accretion rates calculated here (Starrfield et al.|[r985| l. 
These novae outbursts will be observed as super soft X-ray 
sources (kT < IkeV) with luminosities < 10^^ erg s~' and thus 
might be observable in nearby merging galaxies ( |Hernanz"&| 
[Sala 2010). 

Such high mass accretion rates can only be sustained if 
star clusters move slowly with respect to the background gas. 
While there are no stellar clusters observed in the Galactic 
disk which bear these anticipated properties (the relative ve- 
locity of the halo clusters to the interstellar medium is in the 
range of 100 km s"'), observations of several cluster knots in 
the Antennae indicate intracluster relative velocities that com- 
parable to the central velocity dispersions. In these systems, 
it has been shown t hat this merging s tate is, however, short 
lived. At < 20Myr ( [Karl et al.||2010| l. Even if this (slowly 
moving) cold gas phase persists for < 20Myr, we expect at 
least 30% of the white dwarf members to double their mass. 
This increase in mass would in turn lead to a sizable increase 



where e is the efficiency for converting gravitatio nal energy 
into radiation. Extrapolating the Mbh - cr relation (|Tremaine 



et al.| 2002| to the star clusters and dwarf galaxies mod 
eled Dy simulations [A, B, C, D] gives black hole masses 
of Mbh = [0.68,0.68,1.42,1.42] x IO'^Mq and accretion 
luminosities of logL^/Ledd = [10' ", 10-'■'*^ 102■l^ lO-*''^^], 
respectively. While the integrated output from this increased 
mass feeding could, in principle, amount to large factors of 
the black hole's Eddington luminosity, it would probably be 
significantly less because we have, for simplicity, neglected 
feedback. 

We thank Thorsten Naab, Elena D'Onghia, Charlie Con- 
roy, Aaron Romanowsky, Jiirg Diemand, and Lars Hemquist 
for useful discussions. We also thank the referee for his/her 
insightful comments. The software used in this work was in 
part developed by the DOE-supported ASCI/ Alliance Center 
for Astrophysical Thermonuclear Flashes at the University 
of Chicago. Computations were performed on the Pleaides 
UCSC computer cluster. This work is supported by NSF: 
AST-0847563 (J.N. and E.R.), NASA: NNX08AL41G (J.N. 
and D.L.), and The David and Lucile Packard Foundation 
(E.R.). 



in luminosity (e.g. Mestel 19521, which in turn modify the 
white dwarf cooling sequence. 



5.2.2. Accretion onto Central Massive Black Holes 
Most early-type galaxies and spiral bulges are now thought 



to con tain massive black holes in their nuclei (Magorrian et al. 
19981 1. Direct evidence is also now available which supports 



the idea that active galaxies are powered by massive black 
holes. It can then be argued that many if not all galaxies could 
be expected to have undergone an active phase and to pos- 
sess a central massive black hole ( [Haehnelt & Rees|1993j ), al- 
though the presence of intermediate mass black holes in dwarf 
galaxies (and so me globular clusters) r emains controversial 
(lUlvestad et al.l2007; Gebhardt et al. 12005; Baumgardt et al. 

llmai 

Large density enhancements in the nuclei of galaxies, dwarf 
galaxies and star clusters would be accompanied by enhanced 
mass supply to the central black hole. For a typical galac- 
tic halo (model E), this enhanced mass accretion results in a 
significant black hole accretion luminosity: 



(31) 



Fig. 1. — Diagram shows the flow geometiy of non-stationary core poten- 
tials with central mass accumulation. A bow shock is produced by the motion 
of the potential with respect to the ambient gas. The initial asymptotic Mach 
number, sound speed and density of the inflowing gas are denoted here as 
/ioc. Coo and Poo. respectively. Also depicted are the characteristic, r,, and 
the core, fDc, radii, together with a sketch of the expected profiles of density, 
p(r), and sound speed, Cs(r), along the flow axis. 
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TABLE 1 

ASTROPHYSICALLY MOTIVATED SIMULATIONS 



Simulation 


7 
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Mp or M), 


Analytic 


Simulated 


Analytic 


Simulated 




Label 




[kms-'] 




[pc] 


[Mq] 


log(Praax/Poo) 


log(Praax/Poo) 


n [pc] 


rs [pc] 




A" 
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10 


2 


1 
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6.03 


4.32 


2.7 


1.7 


3.8 


B" 


1.01 


10 


4 
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3.12 


1.08 


1.34 


ot 


5.3 




1.01 
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2 
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lO** 
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4.79 


1500 
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4.8 
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100 
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200 


10^ 
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11 


Ot 
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1.5 X 10^ 


1.4 X lO"* 


3.23 


Bf 


1.67 


1580 


2.5 


10*^ 


10'5 


0.5 


0.86 


6x 10^ 


3.6 X 10^ 


4.4 



Young star cluster modeled here by a Pluramcr puleniial. 
^ Globular cluster modeled by a Plummer poienlial. 

Slowly moving dwarf galaxy in cool gas modeled by a Hemquist potential. 

^ Slowly moving dwarf galaxy in hot gas modeled by a Hemquist potential. 

^ Milky Way sized halo moving through hoi gas modeled by a Hemquist potential, 
f 

Simulation of galaxy cluster similar to the Bullet Cluster, modeled by a Hemquist potential. 
^ Those simulations resulted in a detached lailshock. 
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Fig. 2. — Density contours of the flow pattern around point mass, Plummer 
and Hernquist potentials set in motion tlirough an initially uniform medium 
with Poo = 1 and fioc = 2. The point mass potential is implemented here 
with a sink boundary following the method of Ruffert ( 19941 by flooring the 
values of density and pressure to 50% of their background values and set- 
ting the velocities to zero inside the sink radius. Color bars show density cuts 
through thex)'-plane with limiting values: log (pmax/ Poo) = [1.2, 1.1, 1.5] for 
point mass, Plummer and Hernquist potentials, respectively. The equation of 
state used here is adiabatic (7 = 5/3). Also depicted are the characteristic, 
(which reduces to the Bondi radius Rg for a point mass potential), and the 
core, OJc, radii. The Plummer and Hernquist potentials have been constructed 
to yield the same circular velocity peak radius: r\. Note that the difference 
betw een in the Plummer and Hernquist potentials is consistent with equa- 
tion (30l. 
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Fig. 3. — Bow shock patterns produced around a point mass (pink), Plum- 
mer (red) and Hernquist (brown) potentials set in motion tlirough an ini- 
tially uniform medium with /toe = 2.0. Also plotted, for comparison, is 
the bow shock expected for a supersonically moving, non-gravitating object 
(dashed line). In this case the bow shock's opening angle is simply given by 
^'bow = arcsin . The size of the box is 40 times the sink size (used here to 
model the point mass potential), depicted by the black circle. 
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Fig. 4. — Density and radial velocity contours for gas flowing around a point mass, Plummer and Hernquist potential at various relative velocities. We exploit 
the flow symmetries to plot both the density contrast (upper, maroon half) and inward radial velocity profiles (lower, blue half) in the same panel. Color bars 
show density and velocity cuts through the xy-plane in units of log(pniax/ Poo) and p/poo. The maximum densities and radial velocities attained in simulations 
1 through 9 are log(pn,ax/pcx)) = [1.5, 1.2,0.8, 1.4, 1.1, 1.4, 1.7,1.5, 1.1] and (pmax/poo) = [3.9,6.2,6.9,1.9,3.1,6.2, 1.7,2.9,6.4], respectively. The equation of 
state used here is adiabatic (7 = 5/3). 
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Fig. 5. — Density cuts alone the axis parallel to the incoming flow for the 
simulations shown in Figure B] The dotted lines give the analytic solutions 
for a point mass, Plummer anoHernquist potential, while the black, pink and 
red lines show the numerical profiles for /Joo = 1 0. 2.0, and 4.0, respectively. 
The shaded regions depict either the core radii or the sink size. 
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Fig. 6. — Density and radial velocity contours for gas flowing with fiac = 
2.0 around a Hernquist potential in the adiabatic (7 = 5/3) and isothermal 
(7 = 1.01) limits. The maximum densities and radial velocities attained in 
simulations for 7 = 5/3 and 7 = 1.01 are log(pn]ax/Poo) = [1.5,5.1] and 
(/imax//ioo) = [2.9,4.9], respectively. 
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Fig. 7. — Same as Figure |5] but for a near isothermal equation of state 
(7 = 1 .01) Th e dott ed li nes are the analytic profiles derived in ^2] specifically 
equations \23\ and \25\ with 7 = 1.01. The dashed lines are tne isothermal 
form of the density equation, {TTJ. The shaded regions show the size of the 
core. 
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Fig. 8. — Density and radial velocity contours for gas flowing with fioc = 
2.0 around a Hemquist potential with varying in the adiabatic (7 = 5/3) 
limit. Dashed circ les depic t the location of the analytic characteristic radius, 
(equations |19| and (20| ), while solid circles give the approximate location 
of the core radius !I5c . 1 he maximum densities and radial velocities attained in 
simulations 1 thi'ough 6 are log(p„,ax/poo) = [1.8, 1.3, 1.0] and (pmax/Moo) = 
[3.3,3.6,3.2], respectively. 
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Fig. 9. — Expected density enhancements as a function of rs/fDc for a va- 
riety of simulations with c^o = lOkms"' . The triangles (diamonds) are mea- 
sured pmax/poo and rs/Uic from the isothermal (adiabatic) simulations. For 
compariso n, th e lin es a re trends of pmux/ Poo from the analytic results from 
equations \20) and ^23) for Plummer potentials with a variety of masses and 
core radii in Coo = lOkms"' gas. 





{) Inward Radial Velocity ^emaJ^r. 

Fig. 10. — Examples of gas accumulation in star clusters. Density and 
radial velocity contours for models A and B. The maximum densities and 
radial velocities attained in the simulations are log(pn,ax/poo) = [4.32, 1.08] 
and (timax/coo) = [3.8,5.3], respectively. 
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Fig. 1 1. — Examples of gas accumulation in dwarf galaxies. Density and 
radial velocity contours for models C and D. The maximum densities and 
radial velocities attained in the simulations are log(pn,ax/poo) = [4.79,0.03] 
and (Hmax/coo) = [4.8.4.0], respectively. 
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Fig. 12. — Gas accumulation for a Milky Way-like halo moving in a ma- 
jor merging environment. Density and radial velocity contours for model E. 
The maximum densities and radial velocities attained in the simulation are 
log(praax/Poo) = 0.77 and (umax/coo) = 3.23, respectively 
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Fig . 13. — Gas accumulation for a galaxy cluster in a major merger environ- 
ment. Density and radial velocity contours for model F. The maximum densi- 
ties and radial velocities attained in the simulation are log(pmax/Poo) = [0.86] 
and (Wmax/coo) = [4.4], respectively. 
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Fig. 14. — The emission properties of the accumulated gas in a galaxy clus- 
ter (model F) and a young stellar cluster (model A). Emissio n is calculated 
using the cooling curves for an optically thin plasma from Pe res et al.H1982^ 
for chromospheric energy losses in the tempera ture range 4.44 X 10^ < T < 
2 X lO^'X^, the curves of |Rosner et al.|jl978| for coronal temperatures of 
2 X l(f <T < 10* A', and bremsstrahlung emission above 10** K. We assume 
7 = 5/3. That is. the cooling term is not explicitly included in the energy 
equation. The ionization fraction is assumed to be purely coUisional, giving 
Xi = 0.003 for model A and Xi ~ 1 for model F. Model A has a temper- 
ature range in the band kT = 0.6 — 7.2eV, while for model F the range is 
/(:r = 15.0-146.0keV. 
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Fig. 15. — The mass accretion rate fraction, /, of the white dwarf clus- 
ter members for both centrally condensed (solid line) and non-condensed 
(dashed line) compact remnant distributions. The gas density profile for 
model A has been used, with poo = 10^^'gcm"^. To calculate the mass 
accretion distributions we have assumed M„d = IMq and u„a = ISkms"'. 
The red (brown) line was calculating using the numerical density profile for 
fiao = 2.0(0.5). The black line assumes the analytic profile for /loo = 2.0. 



